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We construct a consistency test of General Relativity (GR) on cosmological scales. This test 
enables us to distinguish between the two alternatives to explain the late-time accelerated expansion 
of the universe, that is, dark energy models based on GR and modified gravity models without dark 
energy. We derive the consistency relation in GR which is written only in terms of observables - the 
Hubble parameter, the density perturbations, the peculiar velocities and the lensing potential. The 
breakdown of this consistency relation implies that the Newton constant which governs large-scale 
structure is different from that in the background cosmology, which is a typical feature in modified 
gravity models. We propose a method to perform this test by reconstructing the weak lensing 
spectrum from measured density perturbations and peculiar velocities. This reconstruction relies 
on Poisson's equation in GR to convert the density perturbations to the lensing potential. Hence 
any inconsistency between the reconstructed lensing spectrum and the measured lensing spectrum 
indicates the failure of GR on cosmological scales. The difficulties in performing this test using 
actual observations are discussed. 
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I. INTRODUCTION 

In 1998, astronomers discovered that the expansion 
of the Universe was accelerating, not slowing down as 
expected 0, Q • This late-time acceleration of the Uni- 
verse is surely the most challenging problem in cosmol- 
ogy. Within the framework of General Relativity (GR), 
the acceleration originates from dark energy. The sim- 
plest option is the cosmological constant, first introduced 
by Einstein. However, in order to explain the current 
acceleration of the Universe, the required value of the 
cosmological constant must be incredibly small. Particle 
physics predicts the existence of vacuum energy, but it is 
typically many orders of magnitude larger than the ob- 
served values of the cosmological constant (for a review 
see Ref. 0). 

Alternatively, there could be no dark energy, but a 
large-distance modification to GR may account for the 
late-time acceleration of the universe H! (for a review 
see @, In order to realize the late-time acceleration, 
it is necessary to modify GR at cosmological scales. It 
is really challenging to construct a consistent theory of 
modified gravity (MG) and progress in theoretical physics 
is required in order to have fully consistent models. Even 
if successful theoretical models for MG are constructed, 
we need to distinguish them from dark energy models in 
GR via observations. This is indeed possible if one can 
combine various data sets which measure not only the 
expansion history of the Universe but also the structure 
of the Universe However, it has been argued that by 
fine-tuning the properties of dark energy, it is always pos- 
sible to mimic any MG model and hence it is impossible 
to distinguish between the two possibilities d, [l(| ■ 
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In this paper, we argue that it is indeed possible to 
distinguish between MG models and dark energy models, 
and we propose a consistency test of GR which enables 
us to distinguish between these two possibilities (see also 
Ref. [ll|). From large-scale structure in the Universe, we 
can measure three quantities. One is the density pertur- 
bation which can be measured from the galaxy distribu- 
tion and cluster abundance. The second is the peculiar 
velocity which can be measured from the redshift dis- 
tortions of the galaxy power spectrum and the internal 
dynamics of clusters and galaxies. Finally, we can mea- 
sure the lensing potential from weak lensing of galaxies 
and clusters. We construct a consistency equation in GR 
which can be written only in terms of these observables. 
In this way, it is possible to avoid using specific MG mod- 
els or parametrising MG models in order to test GR on 
cosmological scales. Then we discuss how we can real- 
ize this test using actual observations in the future. The 
most difficult part of the consistency test is that all the 
quantities should be measured at the same time and at 
the same location. However, weak lensing only measures 
the integrated effects of the lensing potential at differ- 
ent redshifts and there is no simple way to reconstruct 
the lensing potential at a given redshift. In this paper, 
we propose a novel way to overcome this problem by re- 
constructing the lensing power spectrum from measured 
matter perturbations at given redshifts. 

Before going into the construction of the consistency 
test, we first summarize the assumptions that we make 
in this paper. First of all, we assume the cosmo- 
logical principle and describe our background universe 
as an isotropic and homogeneous Friedmann-Robertson- 
Walker universe. Thus we do not consider the possibility 
of an inhomogeneous universe which has been extensively 
studied in the literature as a possible explanation for the 
late-time accelerated expansion of the Universe (see [l2| 
for a review). The cosmological principle can be tested 
from observations [H, [13, Q3, Gil and we assume that 
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we can use this to exclude the inhomogeneous universe. 
Next, we assume that the energy- momentum tensor for 
dark matter and baryons is conserved. In MG models, 
we try to explain the late-time accelerated expansion of 
the Universe without introducing any exotic fluids, so 
this is a reasonable assumption. In dark energy mod- 
els, dark matter can be coupled to dark energy (see for 
example However, coupling between baryons and 

dark energy is strongly constrained by local experiments. 
Thus only the conservation of the energy-momentum ten- 
sor for dark matter can be broken. This creates a dif- 
ference between peculiar velocities of dark matter and 
baryons, which can be tested from observations. This 
test will be investigated in a forthcoming paper [l8| and 
we assume that we can use this to exclude the possibility 
of the breakdown of energy-momentum conservation for 
dark matter. The final assumption is that MG models 
are metric theories and there exists a Newtonian limit 
of the theory. This is necessary to reproduce the very 
tight constraints on deviations from GR at solar system 
scales. With this assumption, on subhorizon scales, only 
the time-time component and the anisotropic part of the 
space-space component of the modified Einstein equa- 
tions for perturbations are important. From the Bianchi 
identity and the conservation of the energy-momentum 
tensor, other components of the modified Einstein equa- 
tions should be trivially satisfied. 

The structure of this paper is as follows. In section II, 
we summarize the basic equations for background cos- 
mology and perturbations in GR and MG. In section III, 
the consistency test of GR is constructed, written only 
in terms of observables. We provide a set of consistency 
equations depending on the assumptions for the proper- 
ties of dark energy in GR. In section IV, we propose a 
method to perform the test by reconstructing the weak 
lcnsing power spectrum from measured matter perturba- 
tions. Section V is devoted to conclusions. 



II. BACKGROUND AND PERTURBATIONS 

In GR, the Friedman equation determines the expan- 
sion rate of the Universe 

—^—Pt . Pt —Pm + Pde, 



(1) 



where p m is the matter energy density which include both 
dark matter, and baryon and the dark energy density p de 
accounts for the late-time acceleration of the Universe. 

On the other hand, in MG models, the expansion 
is determined by a modified Friedman equation H 2 = 
F(8irGp m ) where F is a function determined by the 
underlying modified theory of gravity. We can always 
rewrite this modified Friedman equation as 
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For a given MG model, it is always possible to find 
a dark energy model which gives exactly the same ex- 



pansion history, i.e. the same H. This gives the re- 
lation between total energy densities in two models as 
pf*/pF° = G e ff(t)/G. 

We use the Newtonian gauge to describe the metric 
perturbations 



ds 2 



(1 + 2*)di 2 + (1 + 2$)a 2 6 ij dx i dx j . (3) 



The density fluctuation is defined by S — (p — p)/p for 
each component of matter where the quantities with bar 
denotes the background quantity. The divergence of the 
peculiar velocity is — VjTg/(p + P). In GR, the per- 
turbed Einstein equations at sub-horizon scales give 
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(4) 
(5) 



where a is an anisotropic stress. The quantity with the 
subscript t is the total contribution from matter pertur- 
bations 



PtSt = PmSrn + Pde^de, 
(pt+Pt)0t = (pde + Pde)0de + Pn 
(Pt+Pt)<7t = (Pde + Pde)&de- 



(6) 
(7) 
(8) 



As explained in the introduction, we assume that 
there is no interaction between dark energy and dark 
matter/baryons. Then conservation of the energy- 
momentum tensor gives 
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and 
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where SP de is the pressure perturbation of the dark en- 
ergy fluid. 

Assuming a Newtonian limit in the MG model, the 
modified Einstein equations yield the following equations 

El 



fc 2 $ = 4^G e// (a,k)a 2 pf G <5 t MG , 
$ = -r)(a, k)*, 



(12) 
(13) 



where G e //(a, k) and 77(0, k) should be determined by the 
modified theory of gravity. In the MG model, we have 
only dark matter (p t = p m , St = S m ) and the conserva- 
tion equations are given by Eq. ([5]). All other compo- 
nents of modified Einstein equations should be satisfied 
trivially thanks to the Bianchi identity. 
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III. CONSISTENCY TEST OF GR 

First, let us consider simple dark energy models like 
quintessence, where we can safely ignore clustering of 
dark energy, at least under horizon scales. In this case, 
the metric perturbations are given by 



= ^Ga 2 p GR 6g R , 



\$r GR = — $ Gfl 



in the GR dark energy model. In reality there is no 
direct measurement of the curvature perturbation $. 
The weak gravitational lensing provides the combination 
<E>_ = ($ — \f r )/2 which determines the geodesies of pho- 
tons. Then we can define a consistency parameter as 



(o,k) = 



k 2 <f>. 



A-KGa 2 PmSm 



(15) 



where G is the Newton constant measured by local ex- 
periments. If GR holds, c^ 1 ) = 1 (see also Ref. 19] for an 
early attempt to test GR from large-scale structure and 
gravitational lensing). We emphasize that all quantities 
in the consistency relation are observables and we do not 
need to assume any theory. It is easy to see that unless 
(1 + r)- 1 )G eff /2G = 1, we have aW ^ 1 in MG. We 
should emphasize that this test can allow a special class 
of MG with (l+r)- 1 )G eff /2G = 1 to pass the GR consis- 
tency test cr 1 ) = 1. However, if G e ff/G ^ 1, then rj ^ 1, 
which means that the GR relation $ = — *S> is modified. 
In order to check this, we need one more equation which 
relates the Newtonian potential \& to matter-energy per- 
turbations. Since energy-momentum is conserved, the 
continuity equation for matter is not altered and Eq. §§§ 
gives k 2 ^> = d(a8 m )/dt. Then we define a consistency 
parameter for r\ ^ 1 as 



Kk) = -- 



&itGa 2 p m 5, 
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where V m = 2(dd m /dt + H9 m )/{2>aH 2 ). V m can be con- 
structed by direct measurements of the peculiar velocity 
and acceleration of matter. We should emphasize that 
only energy-momentum conservation for dark matter is 
used. Thus even if there is a peculiar velocity of dark en- 
ergy 8de, the reconstruction of the Newtonian potential 
is not affected. If a' 1 ) = 1 and or 1 ' = 1, we can conclude 
that there is no modification of GR and dark energy has 
no clustering. 

However, even if one finds c*W ^ 1 or cr ) ^ 1, this 
does not give definite evidence for the breakdown of GR. 
In general, dark energy can have non-trivial clustering. 
In fact, the conservation equations for dark energy per- 
turbations are not closed and the pressure perturbation 
5Pd e and the anisotropic stress ade should be determined 
by the microphysics of the dark energy fluid. Given that 
we do not know the origin of dark energy, we can treat 
these two functions as arbitrary functions. Then it is 
always possible to mimic $ and $ in a MG model by 
tuning SPde and ade in a dark energy model in GR, be- 
cause Spde mimics G e ff and ade mimics r\ [9(. This fact 



can be clearly seen by rewriting Poisson's equation in GR 
with clustering DE as 



fc 2 3> = AnGa 2 (p m 8 m + p de 5de) 
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(17) 



(14) 

Then, if we only observe dark matter fluctuations 5 n 



cannot distinguish between clustering DE and modified 
gravity as the clustering DE mimics the modification of 
Newton's constant 



G e// (a,k) = G 1 



PdeS, 
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Thus if we find ^ 1 and a' 1 ) 7^ 1, there are still two 
possibilities. One is MG and the other is dark energy 
model with non-trivial clustering. 

However, this argument is clearly based on the assump- 
tion that we only observe the clustering of dark mat- 
ter [ll|. Once the dark energy clusters like dark matter, 
both perturbations participate in forming halos and thus 
galaxies. Thus the galaxy distribution should trace the 
total matter perturbations, which are the sum of the dark 
matter and dark energy perturbations. In the same way, 
the cluster abundance is sensitive to the total matter per- 
turbations [ill ] . When we measure the total matter per- 
turbations from the galaxy distribution, bias between the 
galaxy density fluctuations and the total matter pertur- 
bations becomes the main issue. In this paper, we assume 
that bias can be measured independently by higher order 
statistics or other methods. We will come back to this 
issue in the conclusions. 

Assuming we can measure the total density perturba- 
tions from observations, it becomes possible to improve 
the consistency test and we can distinguish clustering 
dark energy and modified gravity. Including the dark 
energy perturbations, Poisson's equation is now written 
as 



2 JGRsGR 



k ^UH = 47rGa ' p GR S: 

in the GR dark energy model, and 



(19) 



k 2 $ MG = 4irG eff (a, k)a 2 pf G Sf G , (20) 

in MG models. Let us consider two models that give the 
same a^^a, k) ^ 1. In these two models, $ Gfl = $ MG . 
Then the relation between total density perturbations in 
the two models is 
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(21) 



Here G e f / (a) is the effective Newton constant determined 
by the background Friedman equation ([2]) and G e //(a, k) 
is the effective Newton constant determined by the Pois- 
son equation (fT2")) . In general these two effective Newton 
constants are different, as is seen in explicit examples like 
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the DGP [H and f(R) gravity models @. We expect that 
this is a very general feature of MG. 

Then we can define the consistency parameter as 



a 



(2), 



a,k) 



fc 2 $ 



A-KGa 2 p t S t 



(22) 



In GR, a*- 2 -* — 1. The total energy density /?t is recon- 
structed from iJ using the Friedman equation in GR, 
Eq. (fTJ). In this case, we cannot assume the GR relation 
<I> = — iff and we should combine weak lensing that mea- 
sures $_ and peculiar velocity that determines \&. Then 
the consistency parameter is given by 



Ak 2 
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(23) 



Again all the quantities in this equation can be measured 
from observations and we do not need to assume any 
theory to perform the test. From Eq. (|2~Tj) and using the 
fact that in clustering DE and MG, V m and H are 
the same, we find 



,(2) 
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G eff (a,k) 
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(24) 



in MG models. This test (|23|) is essentially the same as 
the consistency relation proposed in [111 ]. Here, we have 
shown that this consistency relation is a test of the differ- 
ence between the Newton constants for the background 
cosmology and for sub-horizon perturbations. 



IV. CONSISTENCY TEST FROM PROJECTED 
POWER SPECTRUM 

All observables in the consistency test should be mea- 
sured at the same time and at the same location. How- 
ever, in reality, there is no simple method to compare the 
statistical observables S t and V m in the consistency 
equation Eq. (|23p . This is because weak lensing is an in- 
tegrated effect of the lensing potential determined by <&_ 
along the line of sight and it is impossible to reconstruct 
the lensing potential at a given redshift. In this paper, 
we propose a novel way to accomplish the test using the 
projected power spectrum (see [20j for a review of the 
projected power spectrum). 

The projected angular power spectra Cf x of any pair 
of perturbations X and X' are given by 



Cf x ' = ^-J dDDW x (D)W x '(D)^A co 



(25) 



where D is the comoving distance and we used the Limber 
approximation which is valid for large I. A^(ao,k) is 
the rms amplitude of curvature fluctuations on comoving 
hypersurfaces at some time given by a = ay during the 
matter-dominated era. The window function W (D) is 
determined by the property of the quantity X as is shown 
below. 



The deflection angle d due to gravitational lensing is 
defined by the gradient field of the lensing potential, d = 
where 



D S ~D 

dD— $_ 

DD„ 



(26) 



and D s denotes the comoving distance to the source 
galaxies distributed on the thin redshift shell labeled by 
s. Then the window function for <t> is 



W' + (D) = -2G$_(a,k) 



(Ds - D) 



(27) 



where the growth function G&_ (a, k) is given by 
G<&_(a.k) — <f>-(a, fc)/$_(ao, k). The angular power 
spectrum of the deflection angle is given by C\ dd = 
l{l + 1)CJ , and Cf dd is related to the convergence 
power spectrum C| KK as Cf KK = t(l + 1)C| dd /A. 

The density perturbations are measured at a given red- 
shift labeled by i at the comoving distance Dt from the 
observer. In the approximation of the quasi-static evolu- 
tion of perturbations, the projected angular power spec- 
trum can be written in discretized form at the given red- 
shift bin i as 

C\ xx ' = ^AAAW^(A)IF X '(A)^A CC . (28) 

The galaxy density fluctuation 5 g is a biased tracer of the 
total density perturbation St , which includes the possible 
presence of dark energy clustering: S g = bSt- The window 
function for the S g component is given by 



r- 



(29) 



where n, is the number density of galaxies at z = Zi, 
rig(zi) and 6j is b(zi). The growth function Q$ t is given 
by 



Gs m (oi, fc) 



a t <I>(a t ,fc) 
$(a ,fc) 



(30) 



if there is no clustering of dark energy (<5 t = S m ) and 



G& t (ai, k) 



» m gg$(q t ,fc) 
a 2 m( ai )^{a ,ky 



(31) 



if there is clustering of dark energy. Both growth func- 
tions are normalized to ao during the matter-dominated 
era. 

The redshift-space power spectrum of galaxies is 
anisotropic due to the peculiar velocities of galaxies. We 
use an approximation that peculiar velocities of galaxies 
trace that of matter m despite the existence of dark en- 



ergy peculiar velocity 9d e [ill] - This allows a statistical 
measurement of the peculiar velocities 9 m (see for exam- 
ple [ll|). It is possible to extract the power spectrum of 
9 m independently of galaxy bias [22]. The derivative of 
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9 m appearing in V m can be derived from a direct mea- 
surment of the acceleration field 23], or estimated from 
measured peculiar velocities at different time slicings [3] • 
Here we assume that it is built from one of these meth- 
ods. Then the window function for V m in each bin i is 



W iV -(D i ) = lg Vm (a i ,k)^-n i 



e 



o u 2 n 2 ' 



where the growth function Qv m (flj > k) is given by 

n m H^(a t ,k) 



Gv m (a,i,k) 



a}H*{ai)y{a ,ky 



(32) 



(33) 



The galaxy power spectra C] SgSg and C\ VmVm can be 
compared directly at each t because these quantities are 
measured at the same time for a given measured comov- 
ing distance Di, and at the same scale k = i/Di in the 
Limber approximation. However, we can observe only 
the integrated effect of the lensing potential and it is im- 
possible to measure at a given time and location. 
Thus we propose a new way to perform the consistency 
test. We statistically reconstruct the lensing power spec- 
trum C| dd by replacing the lensing potential <E>_ with 
the measured density perturbations S t and peculiar ve- 
locities V rn at given redshifts. In the reconstruction, we 
use Poisson's equation in GR to relate 8t to the 3D cur- 
vature perturbations $. Thus if there is inconsistency 
between C\ dd and Cf dd , it indicates a modificaiton of 
the GR Poisson equation, and hence the breakdown of 
GR in structure formation. 

We illustrate how this idea works in the following. The 
continuous contributions from the lensing potential in 
Eq. ([23]) can be discretized 



C; 



dd 



9w 2 - 
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D s D i 



-A*_ s _(ai,fc). (34) 



The first test is to assume that there are no dark en- 
ergy perturbations. Then A$_d>_ (a^, fc) can be replaced 
by the measured matter perturbations from the galaxy 
distribution 



9 Df ( dz 
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(35) 



Substituting this into Eq. (f3"4"|) . we derive the recon- 
structed power spectrum 



i=l 



A ? (A,- A) 2 (dz_ 
Dl \dD' 



-mb, 



-C 



(36) 



In the upper panel of Fig. 1, we demonstrate how ac- 
curately this reconstruction can be done, using LCDM 
with the cosmological parameters tt m = 0.25, = 0.05, 
Hq = 72 km/sec/Mpc, and a flat model with no DE 
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FIG. 1: Upper panel: The precision test between Cg dd (solid 
curves) and C^ dd with a redshift spacing Azi = 0.1 (dash curves) 
for 6 different z s running from 0.25 (bottom) to 2.75 (top). We 
restricted the reconstructed power spectra to scales where the non- 
liner effects give less than 5% changes on the power spectra. Middle 
panel: The first consistency test for GR model (solid curves) and 
MG model -f(R) gravity |25| (long dash curves). Bottom panel: 
The second consistency test for GR model (solid curves), sDGP 
model [H H3 (dash curves), and nDGP model [H Hj] (dotted 
curves). 



clustering. The primordial perturbations are character- 
ized by A c< (a ,k) = S 2 c {a Q ,k n )(k/k n )^ n '^T 2 (k), where 



S((a ,k n ) = 4.52 x 10 



0.95, k n = 0.05Mpc _i , and 



T{k) denotes the transfer function. The bin number of 
source galaxy distributions is running from 1 to 6 between 
z = 0.25 and z — 2.75, with the spacing Az s — 0.5. For 
definitcness, we assume that the galaxy sets come from 
a net galaxy distribution of n g {z) cx z 2 e~^ 2 / 1 ' 5 - ) . The 
biasing b{zi) is fixed to 1. We show the accuracy of this 
reconstruction with diverse binning spacing of Azi = 0.02 
and 0.1. 

We define an estimator corresponding to the test 
Eq. Ull) as 



C s e dd (a^) 



C: 



dd 



(37) 



We can also define an estimator which provides a test of 
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anisotropic stress or 77 = 1 as 



afHia^C, 



iV m V m 



(38) 



If the consistency conditions or 1 ' = 1 and cS 1 ^ = 1 are 
broken, these estimators deviate from 1. Then there are 
two possibilities left: MG or dark energy with cluster- 
ing and/or anisotropic stress. In this case, we should 
reconstruct the lensing potential using 5 g and V m . Then 
A^_ $ _ in the discretized approximation is given by 



f dz . 4 4 



32tt 2 / AD, V dD 



(39) 



The reconstructed deflection power spectra Cf dd {a ty2 " 1 ) is 
given by putting this into Eq. (f3"4")l . Then the estimator 
for the final consistency test Eq. (|23]) is given by 



Kt(a 



(2) 



s dd 



(40) 



In the middle/bottom panel of Fig. 1, we show TZi(a^) 
in a /(i?) gravity model and TZe(a^) in the DGP model 
as examples for MG models. Clearly, the deviation from 
the GR consistency relation is larger than the errors com- 
ing from the reconstruction. Of course, in order to apply 
this test to real observations, there are numbers of prob- 
lems we should solve, such as the precise measurements 
of redshifts, the reconstruction of V m and the measure- 
ments of bias for St. We should also refine our formu- 
lation including redshift distortions. However, it is en- 
couraging that it is in principle possible to distinguish 
between two major possibilities for the late-time acceler- 
ation. It would be crucial to calculate the signal to noise 
ratio in future observations and optimize the detection 
of the deviation from the consistency condition. 



V. SYSTEMATIC ERRORS IN CONSISTENCY 
TEST 

The main challenge in realizing the consistency test is 
to understand bias and the distribution of galaxies. One 
possibility to measure bias is to use higher-order statis- 
tics proposed by Ref. [3(| and applied to the 2dF [3lj . 
Assuming a local bias S g = b\5 m + & 2 £m/2 and gaus- 
sian initial conditions, the measured bias in the 2df is 
b = 1.04 ± 0.11. For the SDSS LRGs, Ref. [H reported 
that the measured bias is b = 1.87 ± 0.07 assuming a 
linear biasing. Ref. [33j gives a forecast for the measure- 
ments of bias from the bispectrum. Assuming primordial 
gaussian fluctuations, a 100 Gpc 3 survey could measure 
61 with a few percent accuracy. We should note that 
there are additional systematic errors coming from the 
assumption on the local biasing. The other method is 



to use the redshift distortions to extract the velocity di- 
vergent power spectrum Pe m e m which we explain below 
[34| . In the first test, we assume that dark energy has no 
perturbations and the growth of structure is determined 
solely by the background cosmology. Then we can recon- 
struct the power spectrum of density perturbations 8 m 
using Eq. ([5]), from which we can measure bias by com 



paring it with the galaxy power spectrum Ps g s g ■ Ref- 34 1 
reported that a full sky survey would measure bias at a 
sub-percent level at z < 2. This includes the systematic 
errors in the measurement of the velocity divergent power 
spectrum discussed below. In general, there are addi- 
tional systematic errors because the assumptions made to 
measure bias in the above arguments should be modified 
depending on the nature of gravity. Thus it is very im- 
portant to understand the effect of MG on bias [36l | . O ne 
possibility is to extend the halo model approach 3?J to 
MG models. Here the modified Newton constant affects 
the critical density S c at which the spherical over-density 
collapses. The study of spherically symmetric collapse 
in modified gravity has just begun [38[ and it is neces- 
sary to understand the relation between G e //(a,k) and 
5 C which eventually determines the halo bias. The un- 
certainty in galaxy redshift distributions also introduces 
inconsistency between the reconstructed lensing power 
spectrum and measured lensing spectrum. This is stud- 
ied in detail in [34j |. It is found that the uncertainty in 
redshift measurements give an extra bias to the measure- 
ments of density perturbations. 

In Ref. [34} , the errors in the reconstructed power spec- 
trum coming from the uncertainty in the measurement of 
bias and redshift distributions are estimated as 



s dd 



where 



Ft = 



AC: 



9 D?(£> s -A) 2 



Fl 2 



Ab, 



1/2 



dz 

Id' 1 



O 2 H 4 
"m-^O rii 99 



(41) 



(42) 



An estimation of errors using the redshift distortion mea- 
surements of bias can be found in Fig. 5 of Ref. [34]. 

In the second test, it is required to measure the power 
spectrum of the velocity divergence. One of the most 
promising ways is to use the redshift distortions. Even 
though the clustering of galaxies in real space to have no 
preferred direction, galaxy maps in redshift space show 
an anisotropic distribution due to the peculiar velocities 
of galaxies. Thus the measurements of anisotropics allow 
us to extract the peculiar velocities. The errors on the 
reconstruction of Pe m e m from redshift distortions are esti- 
mated in Ref. [35[ . Although it is challenging to measure 
Pe m S m from the redshift surveys, it is shown that future 
surveys such as Euclid can give the fractional errors on 
Pe m e m with a few % to 10% accuracy depending on scales 
for k < 0.1 Mpc -1 . This is one of the major systematic 
errors in the second test and we should probably wait for 
the surveys like SKA to achieve sub % accuracy. 
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VI. CONCLUSION 

In this paper, we proposed a consistency test of gen- 
eral relativity by combining geometrical tests and struc- 
ture formation tests using large-scale structure. From 
the geometrical test using supernovae, cosmic microwave 
background and baryon acoustic oscillations, we can re- 
construct the Hubble parameter H. From structure for- 
mation, we measure the density perturbation S t , the pe- 
culiar velocity function V m and the lensing potential 
We constructed a consistency test which is written only 
in terms of these observables, Eq. (|23|) . The major ad- 
vantage of this approach is that we do not need to assume 
any theory for modified gravity models to test GR on cos- 
mological scales. This test is essentially the same as the 
consistency relation derived in [ll| . We have shown that 
this test probes the difference between the Newton con- 
stants in the background cosmology and in large-scale 
structure. The main obstacle in realizing this test was 
that we should measure all the quantities in the consis- 
tency test at the same time and in the same location. 
This is not straightforward because weak lensing mea- 
sures an integrated effect of the lensing potential at dif- 
ferent redshifts and it is impossible to reconstruct <f>_ at a 
given redshift. In this paper, we proposed a way to over- 
come this problem by reconstructing the lensing potential 
from measured density perturbations and peculiar veloc- 
ities. In the reconstruction, we use the Poisson equation 



in GR to reconstruct $ from the density perturbation. 
Thus any inconsistency between the reconstructed lens- 
ing power spectrum and the measured lensing spectrum 
indicates the failure of the Poisson equation in GR. 

In summary, we have shown that, given the as- 
sumptions (i) the cosmological principle, (ii) energy- 
momentum conservation for dark matter and baryon, (iii) 
existence of a Newtonian limit in the MG model, we can 
test general relativity only by using the measured quan- 
tities from the background expansion history and large- 
scale structure of the Universe. We should also emphasize 
that there exist independent ways to test these assump- 
tions from observations. 

We summarized the major systematic errors in the con- 
sistency test in section V. Some of them have been stud- 
ied in detail in Ref. [11], but clearly controlling these 
systematic errors are essential to perform the consistency 
test. 
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